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(57) The invention relates to reservoir evaluation 
and is more specifically directed to a method of 
characterizing the individual response of a layer 
of a multi-layer hydrocarbon reservoir traversed 
by a well, based on down hole flow rate and 
pressure measurements performed during tran- 
sient tests initiated by changes in the surface 
flow rate of the well, the flow rate being 
measured above said layer during one transient 
test and below said layer during another tran- 
sient test. The variations of downhole pressure 
and flow rate with respect to their respective 
values at the initiation of the transient test are 
determined, each of said flow rate variations is 
normalized by the pressure variation aft r the 
sam tim interval within the same transient 
test, thereby to produce a first pressure- 
normalized flow rate function for the level 
above said layer and a second pressure- 
normalized flow rate function for the level be- 
low said lay r, and said first and second press- 
ure-normalized flow rate functions are 
subtractively combined to generate a function 
representative of the individual response of said 
layer. 




Fig. IB Multilayer Transient (MLT) Tot, with SvpitntuJ Presmir* md Flow R*te Mewuie- 
m«mU 
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The invention relates to a method for individually characterizing, from the standpoint of production perform- 
ance, each of the producing layers of a hydrocarbon reservoir traversed by a well. 

An accurate and reliable evaluation of a layered reservoir requires an evaluation on a layer-by-layer basis, 
which involves that relevant parameters, such as permeability, skin factor, and average formation pressure, 
5 can be determined for each individual layer. 

A first conceivable approach for analyzing individual layers is to isolate each layer by setting packers below 
and above the layer, and to perform pressure transient tests, involving the measurement of downhole pressure. 
The layer is characterized by selecting an adequate model, the selection being accomplished using a log-log 
plot of the pressure change vs. time and its derivative, as known in the art. But this method is less than practical 
10 as packers would have to be set and tests conducted successively for each individual layer. 

An alternative approach relies on downhole measurements of pressure and flow rate by means of produc- 
tion logging tools. A proposal for implementing this approach has been to simultaneously measure the flow rate 
above and below the layer of interest, whereby the contribution of the layer to the flow would be computed by 
simply subtracting the flow rate measured below the layer from the flow rate measured above this layer. This 
is in effect would provide a substitute for the isolation of a zone by packers. But this proposal has suffered from 
logistical and calibration difficulties that have thwarted its commercial application. 

A more practical testing technique, called Multilayer Transient (MLT) testing technique, is described by 
Shah et al, "Estimation of the Permeabilities and Skin Factors in Layered Reservoirs with Downhole Rate and 
Pressure Data" in SPE Formation Evaluation (Sept. 1988) pp. 555-566. In this technique, downhole measure- 
20 ments of flow rate are acquired with only one flowmeter displaced from one level to another level. Flow rate 
measurements are thus acquired at different times. However, because fluctuations may occur in the surface 
flow rate, and also because the'change imposed on the surface flow rate to initiate a transient is of arbitrary 
magnitude, it is not possible to determine the contribution of an individual layer by simply subtracting from each 
other the flow rates measured below and above the layer. This complicates the interpretation of test data. 
25 The object of the invention is to enable each layer of a multi-layer reservoir to be characterized on an indi- 

vidual basis from downhole flowrate and pressure transient measurements. 

A further object is to enable such characterization without impractical requirements insofar as acquisition 
of measurement data is concerned being imposed. 

The invention will be made clear from the following description, made with reference to the attached draw- 

30 ings. 

In the drawings : 

- figure 1 A illustrates the isolated zone testing technique, in the case of a three-layer reservoir; 

- figure 1B illustrates the multilayer transient (MLT) testing technique; 

- figure 2 shows an example of a test sequence suitable for evaluating the individual responses of the layers 
35 with the MLT technique; 

- figure 3 is a flow chart describing the method of the invention, with rectangular blocks showing compu- 
tation steps and slanted blocks showing input data for the respective computation steps; 

- figure 4 compares the results of the method of the invention with those obtained from the isolated testing 
technique, based on a simulated example. 

40 In the case of a single-layer hydrocarbon reservoir, well testing techniques allow the properties (permeabi- 

lity, skin factor, average formation pressure, vertical fracture, dual porosity, outer boundaries,... ) of the reser- 
voir - more exactly, of the well-reservoir system - to be determined. A step change is imposed at the surface 
on the flow rate of the well, and pressure is continuously measured in the well. Log-log plots of the pressure 
variations vs. time and of its derivative are used to select a model for the reservoir, and the parameters of the 

45 model are varied to produce a match between modelled and measured data in order to determine the properties 
of the reservoir. 

In the case of a layered reservoir such as the three-layer reservoir shown in figures 1 A and 1 B, a complete 
characterization of the reservoir implies the determination of such parameters as permeability, skin factor, aver- 
age pressure (and others where applicable) for each of the individual layers, because the same model cannot 
so be assumed for all layers. Therefore, such parameters can only be derived from well test data if an adequate 
model can be ascertained for each layer. 

Figure 1 A illustrates the conventional testing technique in which fluid communication between the well and 
the reservoir is restricted to a particular zone isolated by means of packers set above and below this zone, and 
a test is performed by first flowing the well and then shuting it in, and m asuring the variations vs. time of the 
55 pressure in the well during the time the well is shut in. Such a technique allows the response of each individual 
layer to be analyzed, one at a time, since the pressure measured in the isolated portion of the well will only 
depend on the properties of the flowing layer. 

Figure 4 shows simulated pressure and pressure derivative plots vs. elapsed At - the elapsed time for each 

2 



BNSDOCID: <EP 0481866A2J_> 



9 



EP 0 481 866 A2 



10 



isolated zone test starting from the onset of flow. For computing the simulation, the following properties have 
been used for the respective layers : 

Reservoir and Fluid Properties for Simulated Example 

x f (ft) X « r e (£t) 

200 

1.10"* 0.05 200 
50 5-10" 5 0.01 



Layer 


h(ft) 




k(md) 


Skin 


1 


10 


0.20 


300 


3 


2 


15 


0.15 


100 


0 


3 


50 


0.10 


15 





r = 0.4 ft 

w -5 . 

c = 1.10 /psi 

15 t 

with the following definitions : 
h thickness of the layer 

0 porosity 
20 k permeability 

x f vertical fracture half-length 
X interporosity flow parameter 
co storativity ratio 
r e external boundary radius 
25 Figure 4 shows respective pressure and pressure derivative plots for zones 1 , 2 and 3. For instance, layer 

1 is characterized by the pressure and pressure derivative curves in full line. By identifying such features in 
these curves as the slope of the late-time portion, etc, a model can be diagnosed for layer 1. For more infor- 
mation on model selection, reference is made to Ehlig-Economides, C. :"Use of Pressure Derivative in Well 
Test Interpretation" SPE-Formation Evaluation (June 1989) 1280-2. 

so Figure 1 B illustrates an alternative testing technique, called MLT (Multilayer Transient), which makes use 

of downhole measurement of flowrate in addition to pressure. A production logging string, including a pressure 
sensor 10 and a flowmeter 1 1, is lowered into the well. The logging string is suspended from an electrical cabl 
12 which conveys measurement data to a surface equipment, not shown. 

For each test, starting with a change in the surface flow rate, the logging string is positioned above the 
35 layer of interest so that the flow rate measured by the flowmeter includes the contribution from that layer. The 
logging string is kept at this level throughout the test, and is thus caused to operate in a stationary mode. Press- 
ure and flow rate are acquired at a high sampling rate, e.g. every second, during each test. Figure 2 shows 
simulated data illustrating a possible test sequence and the acquired downhole data (with "BHP" standing for 
downhole pressure and "BHF" for downhole flow rate). 
40 A method will now be described whereby a substitute for the single layer responses as obtained by isolated 

zone tests can be derived from MLT test data. 

We assume that transient tests have been performed with the flowmeter respectively above the upper limit 
and below the lower limit of a zone I of the well corresponding to the layer of interest. Evidently, measurements 
acquired with the flowmeter below the lower limit of zone I will also be used as the flow rate measurements 
45 above the upper limit of the zone lying immediately below zone I. 

Let T k . Ti be the start times of the two transient tests, performed with the flowmeter respectively above and 
below the layer of interest, and At the elapsed time within each test. Pressure measurements yield the variation 
of pressure vs. elapsed time : 

APwf(T k + At) for the test starting at T k 
so ApwrfTi + A t ) for the test starting at time T,. 

Flowrate measurements acquired at level J above zone I during the test starting at tim T k yield a flow rate 
variation : 

[Aq(T k + AOL 

Likewise, flow rate measurements acquired at level J+1 below zone I during the test starting at time T, yield 
55 the flow rate variation : 
[Aq(T, + At)] J+1 

We normalize the MLT data obtained during the test starting at T k by forming, for each value of elapsed 
time At|, the ratio of the flow rate variation to the simultaneous pressure variation : 

3 
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PNRj(AW APUT, + At|) 
The same computation yields for the test starting at T| a ratio: 

5 PNR J+1 (5t,)- APwf(Tl + Atl) 

The pressure-normalized ratios pertaining respectively to level J above zone I and level J+1 below zone I 
are subtractively combined to provide a time-dependent data set which characterizes the individual response 
of layer I. 

In the described embodiment, a suitable entity is formed as the reciprocal of the difference between the 
10 ratios PNRj and PNR J+1 : 

RPNR ^~ PNFUAU-PNRj^At.) 
Although the measurements above and below zone I are made at different times and follow changes in 
surface flow rate which may be (and are generally) different in magnitude, the ratios PNRj and PNR^ may be 
15 subtracted because the normalization provides correction for flow rate fluctuations and for the magnitude of 
the flow rate change which has initiated the transient. 

The "reciprocal pressure-normalized rate" (RPNR) pertaining to layer I is a suitable substitute for the press- 
ure change obtained in the context of an isolated zone test. A log-log plot of the RPNR vs. elapsed time thus 
provides a response pattern for the layer of interest. 
20 Likewise, the log-log derivative plot of the RPNR vs. elapsed time provides an equivalent to the pressure 

derivative response obtained in an isolated zone test. 

Superposition effects may have to be taken into account. Superposition effects result from the fact that the 
well has produced at different rates. When the rate is increased from a first value Q1 to a second value Q2, 
the measured pressure drop will be the sum of the pressure change resulting from the change in the rate and 
25 the pressure changes resulting from previous rate changes, including Q1 (see Matthews and Russell, Pressure 
Buildup and Flow Tests in Wells pp. 14-17, Vol. 1 - Henry L. Doherty series, SPE-AIME, 1967). Superposition 
effects may be insignificant if the change in the surface rate is a large increase. However, superposition effects 
may entail gross distortions in the case of a decrease in flowrate, particularly for features pertaining to reservoir 
boundaries. 

30 Correction for superposition involves that derivation of the RPNR be made with respect to a superposition 

time function rather than to elapsed time At In this respect, reference is made to a publication SPE 20550 
"Pressure Desuperposition Technique for Improved Late-Time Transient Diagnosis" by C.A. Ehlig-Economides 
et al. The following description relies upon this work and will refer to the equations presented in this referenc 
as "SPE 20550 Equ." followed by its number. 
35 The RPNR derivative is computed so as to correct for superposition effects, in the manner described below 

in detail with reference to the flow chart of figure 3. 

The result of the computation is the RPNR derivative for every layer. Fig. 4 shows such RPNR derivatives 
for zones 1, 2 and 3 and compares them with the respective single-layer pressure derivative plots which would 
result from the isolated zone test. It is apparent from figure 4 that the RPNR derivative mimics the single-lay r 
40 pressure derivative as regards the meaningful features of the curves (trough, inflection points, line slopes). 

The RPNR and RPNR derivative are thus efficient tools for individually characterizing a given layer i.e. for 
diagnosing a model for this layer. 

It is to be noted that for the RPNR and RPNR derivative to be determined, no specific constraint is impos d 
on the test sequence. The only requirement is that in addition to pressure, measurements of downhole flow 
45 rate variations vs. time are available both above and below the layer under investigation. 

The flow chart of figure 3 provides a detailed description of the steps involved in the computation of the 
RPNR derivative. Rectangular blocks indicate computation steps while slanted blocks indicate data inputting 
steps. 

Input block 20 recalls the above-mentioned definitions of flow rate qj, q j+1 and pressure Pw, measured dow- 
50 nhole during MLT tests. J is the level above the zone of interest, J+1 is the level below that zone. The elaps d 
time variable At, is defined within each transient test, the starting point being the time T k , T,, of change in the 
surfac flow rate. 

Th computations of block 21 provide the pressure change variation and downhole flowrate chang vari- 
ation vs. laps d time. 

55 The respective pressure-normal iz d rates PNR for levels J and J+1 are computed as explained above and 

recall d in block 22. 

Block 23 recalls the computation of th RPNR pertaining to the zone lying between levels J and J+1 , defined 
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as the reciprocal of the difference of the PNR's. 

Input block 24 indicates that the input data for superposition correction (also called desuperposition) are 
the production rate history data : the times of surface rate changes ,.T h the surface flow rates Q(T1 ), Q(T2) .... 
with Q(T1) being the rate from time 0 to T 1( and the downhole flow rates q(T1), etc. 
5 Block 25 gives the expression for the superposition time function t,^, corresponding to SPE 20550 Equa- 

tions (16), (8) brought together. This function is computed for the transient which is considered representative 
i.e. which shows minimal distortion in its late-time period. As explained above, due to superposition, distortion 
will be minimal for the test which starts with the largest increase in surface rate. Block 26 indicates that the 
derivative of pressure variation with respect to the superposition time function tg up is computed for the represen- 
ts) tative transient mentioned above. 

The computation of block 26 yields, for this representative transient, the derivative of pressure change with 
respect to the superposition time function tg Up . From a log-log plot of this pressure derivative vs. elapsed time, 
the slope a of the late-time portion is computed, as indicated by block 27. 

Then, based on the assumption that the pressure change follows a trend represented by 
is Ap^At) = m e (At)a + b 

the slope m e is computed as indicated by block 28 and explained in that portion of SPE20550 which follows 
Equation (21). 

A desuperposition pressure function psup e (At,) is then computed as indicated in block 29, after SPE20550 
Equation (20). 
20 This leads to a corrected pressure change : 

Apwf(At|) - psuPe(Ati) 

Block 30 indicates that the function known in the art as a deconvolution Ap dd , can then be derived from this 
data set. At this point, a choice between two routes must be made depending on the "smoothness" of the decon- 
volution data set Ap dd obtained from the step of block 30. The data will be considered "smooth" if they provide 
25 a definable pattern. If on the contrary, the data are erratic and show no consistent pattern, they are "not smooth". 
Thus block 31 consists of a test as to the "smoothness" of the data set Ap dd (Ati). 

The general expression for the RPNR derivative with respect to In(At) is as follows : 



30 



DRPNRjjCAt.) 
dln(At) 



(RPNR IJ (Ati))' 



dRNP - ( At . ) 2 

± — (PNR,( At . )) 

dln(At) 
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dRNP j+i< A ti> 

dln(At) 



< PNR J + l< At i» 
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If the answer to the test 31 is "Yes", then the RPNR derivative can be computed by substituting the decon- 
volution derivative 

dApdd 
dln(At) 

for the derivative In(At) of the rate normalized pressure RNP(At(), which is the reciprocal to the pressure-nor- 
malized rate PNR. 

This leads to the expression of block 32 for the RPNR derivative. 

If the data are not sufficiently smooth, recourse will be had to the downhole rate-convolved time function 
tsFFtc. expressed by SPE20550 Equ.(24), recalled in block 33. An approximate RPNR derivative can then be 
computed by the expression indicated in block 34, obtained by substituting the corrected convolution deriva- 
tive : 

Psuo(At,)-RNP(A10 

dtsFRC 

for the derivative vs. In(At) of RNP(Atj). 
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Claims 



1. A method of characterizing the individual response of a layer of a multi-layer hydrocarbon reservoir 
traversed by a well, from downhole flow rate and pressure measurements performed during transient tests 

5 initiated by changes in the surface flow rate of the well, comprising the steps of : 

- determining, for each time interval after the initiation of the respective transient, the variations of dow- 
nhole pressure and flow rate with respect to their respective values at the initiation of the transient test, 
the flow rate being measured above said layer during one transient test and below said layer during 
another transient test, 

10 - normalizing each of said flow rate variations by the pressure variation after the same time interval 

within the same transient test thereby to produce a first pressure-normalized flow rate function for the 
level above said layer and a second pressure-normalized flow rate function for the level below said layer, 
and 

- subtractively combining said first and second pressure-normalized flow rate functions to generate a 
is function representative of the individual response of said layer. 

2. The method of claim 1 , wherein the combining step comprises forming the reciprocal of the difference be- 
tween said first and second pressure-normalized flow rate functions. 



20 



The method of claim 1 or claim 2, comprising the step of differentiating with respect to logarithm of time 
said representative function to form a derivative function representative of the individual response of the 
layer. 

4. The method of claim 3, wherein the differentiating step includes correction for the effects of superposition 
25 resulting from changes in the surface flow rate of the well prior to each transient test. 
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Fig. Isolated Zone Test 
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Fig. IB Multilayer Transient (MLT) Test, with Sequential Pressure and Flow Rate Measure- 
ments 
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PNR and RPNR Computehon Tlow Charr 
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Input Att = 0, t\-to,t2~ to," "»'i-to," ' }h- T *,U 
qj(Tk +&U),Pwf(Tk + Ati), V*iCn * P«*f(Ti+&U) For t = 1/ • • , k; ■ ' , i y 



I 
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Compute [Aq(T fc ,£ + Afifl^i = <?J,J + l(7M + Ati) - qj,J+i(Tk,l) 
&pu,f(Tk,L + Aft) = Pmf(Tk,L + Ati) " P"»^ T "iO 
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Compute PNHvlN = AfroKTM^i) 



I 



Compute RPNR difference 



RPNRrrfAtfl- PNRj(At , } . PNRJ , 1( A7T) 




Input- limes for surface rate changes: 

ri,T 2/ - ; T/c, • ,r t; 

and downhole and surface rates : 

qj,J>l(T-,), W+ l(T 2 V • f <KJ«.l(7* I l); Q(Tt),Q(73eV 
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Compute taup for representative transient : 
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ComDute dAp<uF(At ^ For represenl-alive transient 
r dt SU p 




FIG. 3 
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Compute slope, a,of late- time pressure derivahveon plot oF log f^gg^] versus 
log (At). 



Compute slope , m , for same late- h me date on plot of p^f (At) vensu3(A tf for a ± 0, 
or p^At) versus In (At) for a = Q. 



i 



29 



Compute ftup^ l 0 usin 9 E< 1- ( 20 ^ for a *°> or 

^ 
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Compure Hiedeoonvolurion ( Ap,jd3 ( j+i(Af t ), from rhe following systems of equations-. 



32 j- ^lOOTH, 



1 
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IF rhe ApddJ j +1 ( At Jare sufficiently smc»rh,CDmpute rhe RPNR derivative as rot- 

din (AO 



dA pddj^(An) CPNR ^ i(M)) 2j 



33 C 



\ II- the A PddJ,J-»-i(^ r i) are nor sufficienHysmooth, compute t&FRC© • 

34 | 

\_ Gompufe approximate RPNR derivative • 

dRPNRj^A^) (RPNRlj(Ati)) 2 [dlPH UP ,(^) - R " p j(*ft>l (PNPj( Afr)r 
din (At) ^ L dtsFRc 

_ dlp &u0 JAM-RNP J .ifAt«)| (p „ OT i(A , n 2l 

d^SFRCo J 

whereRNPj,,^*) = _J_ v ^ 



FIG. 3 (conh'd) 
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